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Abstract 

, The hadronization cross section in a quark plasma at finite temperature and 

■ density is calculated in the framework of Nambu-Jona-lasinio model with explicit 

chiral symmetry breaking. In apposition to the familiar temperature effect, the 
quark plasma at high density begins to hadronize suddenly. It leads to a sudden 
^r) ' and strong increase of final state pions in relativistic heavy ion collisions which may 

be considered as a clear signature of chiral symmetry restoration. 

OO 

1 Introduction 

o _ 

It is generally believed that there are two QCD phase transitions in hot and dense nuclear 
matter |l|]. One of them is related to the deconfinement process in moving from a hadron 
gas to a quark-gluon plasma, and the other one describes the transition from the chiral 
symmetry breaking phase to the phase in which it is restored. The theoretical evidence 
from the lattice simulations of QCD at finite temperature for these phase transitions 
has prompted experimental efforts to create the new phases of matter in the laboratory 
during the early stages of relativistic heavy ion collisions. Since it is still difficult to 



extract definite information about the phase transitions from the lattice simulations with 
nonzero baryon density, we need QCD models to investigate the phase transitions at finite 
density. From the recent study || based on the instanton model, the phase structure is 
more rich at high density than at high temperature. There exists a new phase of color 
superconductivity when the density is sufficiently high. 

Chiral symmetry plays a crucial role for the behavior of light hadrons, in determining 
not only their static properties in the vacuum, but also the changes therein that are due 
to in-medium effects. One of the models that enables us to see directly how the dynamic 
mechanisms of chiral symmetry breaking and restoration operate is the Nambu-Jona- 
Lasinio (NJL) model applied to quarks[|J. Within this model, one can obtain the hadronic 
mass spectrum and the static properties of mesons remarkably well. In particular, one 
can recover the Goldstone mode, and some important low-energy properties of current 
algebra such as the Goldstone- Treiman and Gellmann-Oakes-Renner relations. Because 
of the contact interaction between quarks that is introduced in the model, there is of 
course no confinement. A further consequence of this feature is that the model is non- 
renormalizable, and it is necessary to introduce a regulator A that serves as a length 
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scale in the problem, and which can be thought of as indicating the onset of asymptotic 
freedom. 

At low baryon density, the phase transition of chiral symmetry restoration in chiral 
limit is of second order, in the case with quarks of nonzero mass it becomes a smooth 
crossover as a function of temperature and density. This picture obtained in the NJL 
model with two flavors is consistent with present lattice simulations and many models||. 
However, at high density several models^], |4[] suggest that the chiral symmetry restora- 
tion transition is of first order. Assuming that this is the case in QCD, this first order 
transition will remain in the case in which explicit chiral symmetry breaking is present. 
This raises the possibility that some consequences of chiral symmetry restoration obtained 
in the chiral limit may be found in the real world with massive pions. Some authors 
have discussed the signatures of a first-order phase transition, such as the event-by-event 
rapidity fluctuations and the associated variations in HBT radii. 

In relativistic heavy ion collisions, nuclear matter is first heated and/or compressed up 
to a region of both a chirally symmetric and a deconfined phase of matter, and then in the 
cooling process the deconfined partons will recombine into the observed hadrons, which are 
mainly pions. It is important to understand how a quark plasma converts into a hadron gas 
at finite temperature and density, and especially what the influence of the chiral symmetry 
restoration transition is on the hadronization process. Although we are aware of the 
disadvantages of the NJL model, the lack of confinement and non-renormalizability of the 
effective 4-fermion interaction, its strengths, the transparent description of the chiral phase 
transition and the binding of quarks into mesons, it is worthwhile to discuss the effect of 
chiral phase transition on the hadronization rate in this model. From the study [0] of the 
hadronization rate for the conversion of a quark-antiquark pair into two pions, qq — > 2tt 
at zero baryon density and in the chiral limit, the hadronization cross section diverges at 
the critical temperature T c of the chiral phase transition. Since the phase transition at 
low density will be fully washed out when explicit chiral symmetry breaking is taken into 
account, the dynamic divergence will no longer exist in the real world. A natural question 
is whether any information about the chiral phase transition at high density is retained in 
the hadronization rate when a nonzero current quark mass is considered. In this paper, 
we investigate at finite temperature and density the hadronization rate for qq — > 2tc in 
the flavor SU (2) NJL model with explicit chiral symmetry breaking. 

The paper is organized as follows. In Section 2, we briefly review the quarks at the 
mean field level and mesons within the random phase approximation (RPA). We discuss 
especially the phase diagrams of the model in temperature-chemical potential (T— /i) plane 
and temperature-baryon density (T — n b ) plane. In Section 3, we study the behavior of 
the hadronization cross section to first order in a 1/N C expansion as a function of the cm. 
energy s, temperature T and baryon density n^. We focus our attention on the difference 
between temperature and density effects, and the influence of the chiral phase transition 
on the hadronization rate. We summarize in Section 4. 
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2 Phase Diagrams 



The two-flavor version of the NJL model is defined through the Lagrangian density, 

L NJL = - m )V + G{(^) 2 + ferV) 2 ], (1) 

where only the scalar and pseudoscalar interactions corresponding to a and ir mesons, 
respectively, are considered, ip and ip are the quark fields, r is the SU(2) isospin generator, 
G is the coupling constant with dimension GeV~ 2 , and mo is the current quark mass. 

The key quantity for dealing with the thermodynamics of a variable number of particles 
is the thermodynamic potential Q. In the mean field approximation, which is the lowest 
order 0(1) in the 1/N C expansion, it is|| 

nfrr ^ (m g - m ) 2 f d 3 p 

tl(T, fi, m q ) = — 2iYyV> J j^E p 

-2N c N f T / — ^ In 1 + e~^" 1 + e~^" , (2) 



(2tt)* 

where E 2 = p 2 + m 2 , /x is the chemical potential, and N c , Nf are the number of colors 
and flavors respectively. Equation (Q) holds for any value of the quark mass m q , which is 
the order parameter of the chiral phase transition. In the spirit of thermodynamics, the 
physical system is described by that value of m q for which Q(T, /i, m q ) is a minimum, 

80 ^ #Si >0. (3) 



dm q dm q 

The above equality is the socalled gap equation to determine the constituent quark mass 
as a function of T and fi. The inequality is especially necessary in high density region 
where the thermodynamic potential has two minima and one maximum. 

In the chiral limit when the current quark mass = 0, the solution of the gap 
equation separates the region of chiral symmetry breaking with m q ^ from the region of 
chiral symmetry restoration with m q = 0. The phase-transition lines demarcating these 
two regions are given in Fig. la in the T — fi plane and Fig. lb in the T — n& plane. The 
relation between the baryon density n& and the chemical potential fi is rather complicated 
and is also a function of T, 



N c N f f d 3 p 
n b = — — ' 



(27T) 



tanh 7^{E P + fi) - tanh ^{E p - fi) 



(4) 



where the factor 3 reflects the fact that three quarks make a baryon. 

The order of the chiral phase transition is determined by the behavior of the order 
parameter across the phase-transition line. If the quark mass m q is continuous, the phase 
transition is of second order indicated by dashed lines in Fig.l, otherwise it is of first order 
shown by thin solid lines. The tricritical point P which separates the first- and second- 
order phase transitions is located at T = 0.079 GeV and fi = 0.28 GeV or n b /n = 1.95 
with n = 0.17/ fm 3 being the baryon density in normal nuclear matter. Its positions in 
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the T — fi and T — n b planes are close to that calculated in the Landau 6 model and the 
instanton model |J. 

At any critical point on the first-order phase-transition line in the T — [i plane, the 
thermodynamic potential as a function of the quark mass has two minima and one max- 
imum between them. The coexistence of two minima indicates a mixed phase of chiral 
symmetry breaking and chiral symmetry restoration. The two spinodal points on the left 
hand and right hand sides of the maximum are characterized by a stability equal to zero, 



Between the spinodal points, J-^- < 0, the system is unstable with respect to a small 

fluctuation in quark mass. Outside the spinodal points and inside the minima, J^p- > 0, 
the system is metastable. In these regions a finite fluctuation in m q may be amplified, 
and the system may not return to its original state after a sufficient deviation. Outside 
of the minima, > 0, the system is stable. The unstable and metastable regions are 
marked by US and MS in Fig. lb. 

In the real world with nonzero current quark mass, the second-order phase transition 
becomes a smooth crossover and the tricritical point P becomes a critical end-point E of 
a first-order phase-transition line. The end-point E is a second-order critical point and 
shifted with respect to the tricritical point P towards larger fi in the T — ji plane and 
towards smaller rib in the T — rib plane. The first-order phase-transition lines are denoted 
by thick solid lines in the phase diagrams. Their shift with respect to the case of chiral 
limit is due to the fact that it is more difficult for a phase transition to happen in a system 
of massive particles than that in the corresponding system of massless particles. 

In the calculations of the phase-transition lines shown in Fig.l, the 3 parameters in 
the NJL model, namely the coupling constant G, the momentum cutoff A and the current 
quark mass mo, are fixed by the constituent quark mass m q = 0.32 GeV, the pion decay 
constant f n = 0.093 GeV and the pion mass m* = 0.134 GeV (in the chiral limit m^ = 0) 
in the vacuum (T = fi = 0). In the self-consistent mean field approximation, it is well 
known that the masses of the 7r and a mesons satisfy the dispersion relation @] 

l-2GU M (T,fj,,m M ) = , (6) 

where uim is the mass of the respective meson and Um is the associated retarded polar- 
ization function, 



r d 3 p 1 E 2 l l 

II* (T, /i, m n ) = N c Nf J — -— _ ? 4 ^tanh ^(E p + fj) + tanh ^(E p - y) 



{2^fE p El-ml/A 



1 ,„ , , 1 



r d 3 p 1 E 2 — m 2 r 1 1 
II CT (T, fi, m a ) = N c Nf J E 2 -m 2 lA r h 2T ^ p + ^ + tanh ~ ^ 



{2-KfE p E 2 p -m 2 jA 

The pion coupling constant g w is related to the residue at the pole of the pion propagator, 



_ 2 _ dnj^j^m^) 
' " dm 2 
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N c N f r d 3 p E v r 1 ,„ . , 1 .„ ,1 

= ~r / - ^ r 2r (£ ' + ri + ta " h 2f ( ^ - ">] • (8) 

It is easy to see that when > 2m g , ^ = 0, the pions will decay into quark-antiquark 
pairs, 7r — > gg. Therefore, we can determine a hadronization line by 

m^iTh^h) = 2m q (T h ,fi h ) . (9) 

For temperature T > Th and chemical potential /i > /j/, or baryon density rib > n^h, 
mesons become resonant states and no more pions exist in the model, and one has a 
system of interacting quarks. For T <Th and \i < fih or rij, < nbh, a quark- meson plasma 
is formed. This means that in the cooling process of relativistic heavy ion collisions the 
quark plasma begins to hadronize at the point (T h ,n bh ). The coexistence of quarks and 
mesons at low temperature and density is of course an artifact of the NJL model. 

In the chiral limit, the equation @ for determining the hadronization line is fully 
equivalent to the equation for defining the chiral phase-transition line. Therefore, the two 
lines coincide in the phase diagrams. When explicit chiral symmetry breaking is taken into 
account, the two equations become different. While there is no longer a phase transition 
at high temperature, one can still define the hadronization lines via equations (|9p and 
which are shown in Fig.l. 



3 Hadronization Cross Section 

The total hadronization cross section of a single quark in a flavor 577(2) model is defined 
as0] 

this is a measure of how quickly a quark hadronizes from a charge symmetric plasma. 
Since the NJL model is a strong-coupling theory, perturbation theory is inapplicable, and 
we require a selection procedure for the relevant Feynman graphs. We choose the 1/N C 
expansion, and work to the first order. To 0(l/n c ), the relevant T-matrix amplitudes 
for the reaction qq — ► 2n are sketched in Fig.2. It depicts the s-channel amplitude T s 
in which an intermediate cx-meson is produced, and the t-channel which displays a quark 
exchange. It is understood that the crossed diagrams are included for a process in which 
the final state particles are identical, e.g. for uu — > 2n . 

Using an obvious notation for the T-matrix amplitudes, one has for the process uu — > 
2tcq the s-channel amplitude 

2G 

T s = u(p 2 )«(pi)- , ■ rglAwripx +p 2 ,P3)S clC2 , (11) 

1 - 2GU a (px +p 2 ) 

where u and v are quark and antiquark spinors, the symbol S C1C2 refers to the color degree 
of freedom, and A a7cn represents the amplitude of the triangle vertex a — > 7r7r. By assuming 
the qq system to be at rest in the heat bath and evaluating the Matsubara sum in the 
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triangle, one obtains the explicit dependence of the amplitude A a7T7T on the Mandelstam 
Lorentz scalar s = (pi +P2) 2 , and the temperature T and chemical potential /i, 

A m(s , T, „) = W/ j ^M^Jft^l x 

8(p ■ p 3 ) 2 - (2s + 4m 2 )p ■ p 3 + s 2 /2 - 2sE 2 
(s - AE 2 ) {{ml - 2p • p 3 ) 2 - E 2 s) 

where p 2 = \J ' s/A — m 2 , and the Fermi-Dirac distribution function /f(x) — (l + e x ^ T ^j 
In terms of the other two Mandelstam variables t = (pi — p 3 ) 2 and u = (pi — Pi) 2 with the 
constraint s + t + u = 2m 2 + 2m 2 , the direct and exchange contributions to the t-channel 
amplitude are 

T ' = v(p 2 )75 — — \ l5u{pi)gJ Cl c 2 , 

t m q 

T t,exc = ^(p 2 ) 75 7M(Pl " PA \ + mq Ku( Pl )g 2 J cl c 2 ■ (13) 

U — Tflq 

Since the s-channel crossed graph equals the direct term, the differential cross section can 
be written as 

,/ """- 2T0 (s, t, T, //) = - - 1 - 2 ^ £ \2T s (s, T, + 7*»{t, T, /i) + T^(u, T, /i) | 2 , 



dt ' IQits(s — Am 2 ) 



q) c,s 



(14) 

where the prime on the summation indicates the average over spin and color degrees of 
freedom of the quarks in the initial state. 

It is necessary to note that in the s-channel the triangle amplitude A CT7r7r and the a 
polarization function 11^ are both complex functions. Their real and imaginary parts are 
separated from each other by using the decomposition = + m5(x) for A an7T (s + 
ie, T, fi) and H a (s + ie, T, fi). 

After performing the integration over the momentum transfer t, the integrated cross 
section is obtained from 

rtmax d(J - o / \/~S \ ^ 

(T m -,2^ (s,T,fi) = jf dt ^ no (s,t,T,p,) \1 + f B (—)j , (15) 



where we have the limits t max = — |(s — Am 2 ) — |(s — 4m 2 ) + ~y (s — 4m^)(4 — m 2 ), and 
tmin = -\{s - Am 2 ) - \{s - 4m 2 ) for the process uu -> 2n and t min = -|(s - 4m 2 ) - 



\{s — Am 2 ) — |y ( s — Am 2 ) (A — m 2 ) for the other two processes. In equation (|I~5|) we have 

taken into account the Bose-Einstein distribution function fs(x) = (e x l T — l) for the 
final state pions. 

In the same way as equation ( |T4"D was constructed, we can write the differential cross 
sections for the processes uu — > 7r + 7r_ and ud — > 7r + 7r explicitly as 

''""""^ ts,t,T,») = J— 2 ^ \T s (s,T,^) + 2T t ^(t,T,^)\ 2 , 



' UU — >TXJ r 'K 

dt ' 16tts(s — Am 2 ) 



qj c,: 
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% ± ^^") - ^(L^ t^T^T,,)? . (16) 

The total hadronization cross section discussed above in the whole T — // (or T — n^) 
plane and with explicit chiral symmetry breaking was studied in the chiral limit and 
along the T axis in Ref. [^]. The main result is the singularity structure of the cross 
section: a q (s,T) displays singularities at the threshold energy s t h = 4m 2 for any T and 
at the critical temperature T c for any s. The singularity at sth results from the fact 
that the entrance channel has a threshold energy sth, while the exit channel has s = 0. 
The singularity at T c arises from the chiral symmetry restoration reflected on the quark 
exchange in the t-channel: As T — > T c , m q — > 0, and then the amplitude T t,d%r — > oo. 

What we want to stress in this paper is the baryon density effect and the m n de- 
pendence of the hadronization rate. Fig. 3 shows the energy dependence of the inte- 
grated total cross section as a function of energy s for three values of the density n b 
at zero temperature (Fig. 3a) and for three values of the temperature T at zero density 
(Fig.3b). Because of the three-momentum cutoff A, the energy range is restricted to 
4 (Max(m q , m n )) 2 < s < 4 (^(Max(m q , m^)) 2 + A 2 ). At low densities, there is still a kine- 
matical divergence at s t h = 4m 2 , when m q > m n . For ^ = 1.1 shown by thin solid line in 
Fig. 3a, the cross section drops down from the infinity at s t h = 1.33 GeV 2 and then goes 
up again and reaches the second infinity at s = m 2 which is the contribution from the a- 
exchange in the s-channel. Then the cross section decreases and is fairly constant for large 
energies. With the increase of the density, when m n > m q , the threshold energy of the 
reaction is larger than the entrance energy, and therefore the first singularity disappears. 
In Fig. 3a this is the case for ^ = 1.25 denoted by dashed line. Further when the density 
is sufficiently high that the threshold energy s = 4m 2 of the reaction is larger than m 2 , 
the second singularity disappears too, the cross section has no further singularities, see 
the thick solid line for ^ = 1.5 in Fig. 3a. The energy dependence of the cross section at 
finite temperature is similar and shown in Fig. 3b. 

The energy dependence of the singularity structure discussed here is very different 
from that in the chiral limit where there is always only one singularity at the threshold 
s th = 4m 2 . 

The temperature and density effects on the total cross section for three values of 
energy s are given in Fig. 4. At low density, the chiral phase transition is washed out in 
the case of a non- vanishing current quark mass mo = 0.005GeV , and so is the singularity 
in the cross section at the transition point. Even an expected maximum due to the small 
quark mass at this point does not exist. It reflects the chiral phenomenon that when 
the quark mass drops down with the increase of the temperature, the pion mass goes up 
simultaneously, and therefore it becomes more and more difficult for the quark plasma to 
hadronize into pions. At low temperature, while the first-order phase-transition remains, 
the dynamical singularity in the cross section is still washed out since the finite quark and 
pion masses. However, very different from the temperature dependence shown in Fig. 4a 
where the quark plasma hadronizes smoothly and reaches the maximum rate within a 
temperature range AT ~ 0.1 GeV, there is a cliffy boundary in the T — rib plane (Fig.4b) 
separating the quark-meson plasma at low density from the quark plasma at high density. 
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This feature arises from the discontinuities of the quark and pion masses at the first-order 
phase-transition point. When the density decreases, the quark plasma first meets the 
hadronization point and begins to hadronize, and then it reaches the first-order phase 
transition, the sudden increase of the quark mass combined with the sudden decrease of 
the pion mass at the transition point leads to a jump in the cross section. 

At zero temperature, the hadronization point and the chiral transition point are very 
close to each other, see in Fig.l, the quark plasma reaches the transition point immediately 
after the beginning of hadronization. Therefore, we see in Fig. 4b a precipitous jump 
from zero to a maximum for any s. The maximum, although it is not strong, is also a 
consequence of the discontinuities of the quark and pion masses at the transition point. 

4 Conclusions 

The present investigation of the hadronization of a quark plasma in the NJL model to the 
first order in 1/N C has been motivated by two aspects: 1) The chiral symmetry restoration 
in the chiral limit and its consequence in hadronization rate, a dynamical singularity at 
the critical point, are fully washed out when a no n- vanishing current quark mass is taken 
into account. An important question is then what the influence of chiral dynamics is in 
the hadronization in the real quark world. 2) While the chiral phase transition happens at 
both high temperature and at high density, the order of the phase transition is different in 
two cases. The second-order phase transition due to temperature effect becomes a smooth 
crossover, but the first-order transition due to density effect is retained when the explicit 
chiral symmetry is considered. If such a first-order transition occurs in relativistic heavy 
ion collisions, we need to know what its experimental signatures are. 

The NJL model may be used as a model of QCD as long as one is interested in 
the chiral properties of the quark world. In this model the hadronization point of the 
interacting quarks is above the chiral phase transition point in the real world at high 
density, and the two points approach each other when the current quark mass goes to 
zero. Above the hadronization point one has a deconfined phase of quarks, but a quark- 
meson plasma below the point. Therefore, in the process of hadronization of the quark 
plasma, we expect that the retained first-order chiral transition at high density may still 
play an important role in the hadronization rate. 

We have calculated the energy, temperature and density dependence of the total 
hadronization cross section of a single quark in the framework of flavor SU(2) NJL 
model. Corresponding to the motivations mentioned above, our main conclusions are: 
1) The hadronization cross section is normally strongest for small values of the energy 
of the qq pair, and the singularity structure at low energies depends on the temperature 
and density. At low temperature or density, there are two singularities at the threshold 
energy and at the pole position (ra£) of the a propagator. As the temperature or density 
increases, the kinematical singularity at the threshold energy vanishes first when the pion 
mass m n exceeds the quark mass m q , and then the singularity at the pole disappears too 
when the threshold energy of the hadronization process, s = 4:m^, is larger than vp? c . As 
for the critical scattering, the divergence at the critical point for all values of s in the chiral 
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limit is fully washed away. 2) Very different from the temperature effect which leads to 
a gradual hadronization during the cooling process of the quark plasma, the hadroniza- 
tion cross section is fairly constant in the process of hadronization at low temperature. 
The sudden jump of the cross section at the beginning of the hadronization means that 
the first-order transition of chiral symmetry restoration at high density is still reflected 
in the hadronization in the real quark world. A direct consequence of this feature is a 
suddenly strong increase of pion number in the expansion process of relativistic heavy ion 
collisions, if chiral symmetry is restored in the collisions. Our qualitative conclusions here 
derived from the NJL model may be more general because the first-order transition of 
chiral symmetry restoration looks like an essential feature of a strong interacting system 
at high density. 
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Figure Captions 

Fig.l: The phase-transition line and hadronization line in T — \x plane (Fig. la) and 
T — rib plane (Fig. lb). The dashed lines indicate a second-order transition while the thin 
solid lines indicate one of first order in the chiral limit. The point P is tricritical. The 
thick solid lines with critical end-point E are the lines of first-order transition with non- 
vanishing current quark mass. For the first-order transition with and without current 
quark mass, the mixed region is divided into regions of instability (U S) and metastability 
(MS). The hadronization line separating the region of quark plasma from that of quark- 
meson plasma is indicated by thin solid lines in the two planes. 

Fig. 2: The Feynman diagrams for the hadronization of a qq pair into two pions to the 
first order in 1/N C expansion. The solid lines denote quarks and antiquarks, the dashed 
lines denote mesons. 

Fig. 3: The total hadronization cross section as a function of the collision energy s, 
for three different baryon densities (Fig. 3a) at T = and for three different temperatures 
(Fig.3b) at n b = 0. 

Fig. 4: The total hadronization cross section as a function of the temperature (Fig.4a) 
at rib = and as a function of the baryon density (Fig.4b) at T = 0. In each case, the 
collision energy takes three different values. 
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